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We demonstrate for the first the existence of electrically charged BPS vortices in a Maxwell-Higgs 
model supplemented with a parity-odd Lorentz-violating (LV) structure belonging to the CPT-even 
gauge sector of the standard model extension and a fourth order potential (in the absence of the 
Chern-Simons term). The modified first order BPS equations provide charged vortex configurations 
endowed with some interesting features: localized and controllable spatial thickness, integer flux 
quantization, electric field inversion and localized magnetic flux reversion. This model could possibly 
be applied on condensed matter systems which support charged vortices carrying integer quantized 
magnetic flux, endowed with localized flipping of the magnetic flux. 
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Since the seminal works by Abrikosov |l| and Nielscn- 
Olesen Q showing the existence of electrically neutral 
vortices in type-II superconducting systems and in field 
theory models, respectively, these nonperturbative solu- 
tions have been a theoretical issue of enduring interest. In 
the beginning 90's, vortices solutions were studied in the 
context of planar theories including the Chern-Simons 
term §| which turned possible the attaining of electri- 
cally charged vortices |J] also supporting BPS (Bogo- 
mol'nyi, Prasad, Sommcrfcld) solutions [5j, related with 
the physics of anyons and the fractional quantum Hall 
effect Q. In addition, charged BPS vortices also were 
found in the Maxwell-Chern-Simons model @- Fur- 
ther studies were performed involvin g no nminimal cou- 
pling [§[ and new developments ffjl. Il0||. Generalized 
Chern-Simons vortex solutions were recently examined 
in the presence of noncanonical kinetic terms (high or- 
der derivative terms) II ill usually defined in the con- 
text of k- field theories fl2|-ll4i . These k-defects present 
a compact-like support [151 ] . an issue of great interest 
currently Lately, in the context of effective field 

theories, it has also been demonstrated the existence of 
charged BPS vortices in a generalized Maxwell- Chern- 
Simons- Higgs model fl7| . 

Lorentz-violating (LV) theories have been under atten- 
tion in the latest years. The general theoretical frame- 
work for studying Lorentz-violation effects is the stan- 
dard model extension (SME) [l8| which encompasses 
Lorentz-violating terms in all sectors of the minimal stan- 
dard model. In particular, the abelian gauge sector of 
this model is composed of two sectors, a CPT-odd and a 
CPT-even one. The CPT-even part is described by nine- 
teen parameters enclosed in the rank 4 tensor, (Kp) >J ' au ^ , 
endowed with a double null trace and the symmetries of 
the Riemann tensor fl9j . being investigated in many re- 
spects [20j . Lorentz violation was also considered in con- 
nection with the formation of topological defec ts l2ll |22| , 
with particular interest in the Higgs sector [23|]. Re- 
cently, it has been investigated the formation of stable 



uncharged vortices in the context of the nonbircfringent 
Lorentz-violating and CPT-even Maxwell-Higgs electro- 
dynamics, also including LV terms in the Abelian Higgs 
sector of the SME, with new interesting results [24j . 

Up the moment it is known that abelian charged vor- 
tices are only defined in models endowed with the Chern- 
Simons term. This remains valid even in the context of 
highly nonlinear models, such as the Born-Infield elec- 
trodynamics [25| . In this letter we report for the first 
time the existence of abelian charged BPS vortices in a 
Maxwell-Higgs electrodynamics deprived of the Chern- 
Simons term and endowed with CPT-even LV terms. 
This achievement is ascribed to the CPT-even electro- 
dynamics of the SME, whose parity-odd coefficients en- 
twine the electric and magnetic sectors [III [27j in analogy 
to what happens in the models containing the Chern- 
Simons term. The BPS solutions are attained by consid- 
ering a particular fourth-order potential, and can be in- 
terpreted as vortex solutions in a dielectric medium [28[ . 
The charged vortex solutions are localized, having spatial 
thickness controlled by the LV parameter, and present- 
ing localized magnetic flux and electric field reversion in 
some radial region. This phenomenon could be of interest 
in condensed matter physics, mainly in connection with 
superconductivity, particularly with two-components su- 
perconducting systems [29 ]. 



I. THE THEORETICAL FRAMEWORK 



The theoretical environment in which our investiga- 
tions are developed is a modified Maxwell-Higgs model 
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defined by the following Lagrangian density 



II. BPS CONSTRUCTION 



A+3 = -\f^ v - \k» v F w F v " + \D^\ 2 

+ ^(i + K 00 )d^d^ + ^d^d v ^ (i) 

containing a convenient potential U (\4>\ The two 

first terms in ([T]) define the nonbircfringent and CPT- 
even electrodynamics of the SME, whose nine LV non- 
birefringent parameters are enclosed in the symmetric 
and traceless tensor, , defined as 



= (KfY 



(2) 



where (Kp)^ av ^ is the CPT-even gauge sector of the 
SME, whose properties were much investigated since 
2001 [HHI. The Higgs field, 0, is coupled to the gauge 
sector by the covariant derivative, D^cj) = d^<f> — ieA^cj). 
The neutral scalar field, ^, plays the role of an auxiliary 
scalar field, and is analogous to the one that appears in 
the planar Maxwell-Chern-Simons models endowed with 
charged BPS vortex configurations @,[loj]. Note that the 
Lorentz- violating tensor, k^ v , also modifies the kinetic 
term of the neutral field. The potential in Eq. §3§ which 
assures the attainment of BPS first order equations is 
defined by 



V i 



ev 



2(1-8) 



(3) 



observe that it is a fourth-order one. Here, v is the 
vacuum expectation value of the Higgs field, while s = 
tr (Kij). The potential U possesses a nonsymmetric min- 
imum, = 0, \4>\ = v, which is responsible by pro- 
viding topological charged vortex configurations. We 
now analyze static solutions for the model projected 
in the xy— plane, where it recovers the structure of a 
(1 + 2)— dimensional theory. Remembering that in 
(1 + 2)-D it holds F t j = e tj B, F ai = E l (the magnetic 
field becomes a scalar), we write the stationary Gauss's 
law of the planar model related to ([!]) by 



LijdidjAo + eijKQidjB — 2e 2 Ao 



(4) 



In this section, we focus our attention on the de- 
velopment of a BPS framework Q consistent with the 
second order differential equations obtained from the 
(1 + 2)— dimensional version of the Lagrangian ((T|). We 
begin writing the energy density E in the stationary 
regime as 

£ = \{l-s)B 2 + U{\<P\^) + \D^\ 2 
+ ±Lii (diA ) (djAo) + \u, (d^) (dj*) (6) 

+ e 2 A 2 \^\ 2 + e 2 ^ 2 \^\ 2 . 

In order to attain the first order differential equations, 
we first impose the following condition on the neutral 
field 



* = tA , 



(7) 



which is similar to the ones appearing in the BPS vortex 
configurations of the Maxwell-Chern-Simons model 0, 
0J3. By substituting © in Eq. @, we attain 



E = -(l-s)B 2 



e \cf>\ 2 ± eijKoidjA 
2(1-*) 



(8) 



\Dj<t>\ 2 + Lij (cUo) (djAo) + 2e 2 A 2 \</>\' 



After converting the two first terms in quadratic form 
and by using the identity, 

{D^l 2 = \D ± cp\ 2 ± e \4>\ 2 B ± ^e ab d a J b , (9) 

where J b is the spatial component of the current J M = 
i [(j) (D^cj))* — (jfD^(j)\ 1 the energy density takes on the 
form 



£ =2 (1 - S) 



ev 



e |0| ±6ijKoidjA 



(1-8) 



\D ± 0\ 2 ±ev 2 B±±e ab d a J b 



(10) 



where 



= (1 + k 00 ) Si 



(5) 



carries on the CPT-even and parity-even LV parameters. 
We note that is the parity-odd parameter Koi that couples 
the electric and magnetic sectors (26l . [27j , making possi- 
ble the existence of charged vortex configurations which 
are attained even in absence of the Chern-Simons term. 
Otherwise, for K 0i = 0, the temporal gauge A = solves 
the Gauss law, yielding compactlikc uncharged vortex so- 
lutions in. 



+ L^ (diA ) (djA ) + B d J A + 2e 2 A 2 \ ( f ) \ 2 . 

Now, we use the Gauss's law to transform the last three 
terms in a total derivative, rewriting the energy density 
as 



ev 2 - e \<f>\ 2 ± eijK 0i djAo 



(1-8) 



(11) 



\D±</>\ 2 ±ev 2 B + d a Ja, 
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where 

Ja = i^abJb + L a bA dbA + f-baKobB A . (12) 

This energy density is minimized by requiring that 
the squared terms be null, establishing the two BPS con- 
ditions of the model: 



\D±</>\ = , 



B 



± ( eif - e \4>\ j + eijKoidjAo 



(13) 



(14) 



Under these BPS conditions, we attain the BPS energy 
density, 



Ebps = ±ev 2 B + d a Ja, 
implying the total BPS energy, 

E B ps = ±ev 2 / d 2 rB = ±ev 2 ^, 



(15) 



(16) 



which is proportional to the magnetic flux. It is worth- 
while to note that the second term in (|T5|) does not con- 
tribute to the total BPS energy, once the fields go to zero 
at infinity. 



III. CHARGED VORTEX CONFIGURATIONS 

Specifically, we look for radially symmetric solutions 
using the standard static vortex Ansatz 

Ae = - Q(r) ~ n , <t> = vg (r) e m \ A = w(r), (17) 
er 

that allows to write the magnetic field as 

(18) 



a 
er 



The scalar functions a(r) , g (r) and u (r) are regular at 
r = and r = oo, satisfying the appropriate boundary 
conditions: 

g (0) = 0, a (0) = n, J (0) = cte, (19) 
g(oo) = 1, a(oo) = 0, w (oo) = 0, (20) 

with n being the winding number of the vortex solution. 
The boundary conditions above will be demonstrated ex- 
plicitly in the remain of the manuscript. 

We now introduce the dimensionless variable t = evr 
and implement the following changes: 



9 (r) -> 9 (t) , a(r)^a(t),u (r) w (t) 
B^ev 2 B{t), £ -^v 2 £ (t) . 



(21) 



Thereby, the BPS equations (|13I14[) and the Gauss's law 
Q are rewritten in a dimensionless form 



-/ , a 9 
' y t 

5= _0f_ = ± (1 - g 2 ) - Hjjj' 

t (1-s) 



(22) 
(23) 



(tCj')' (tB)' 
(l + \ r )^^--K^--2g 2 Cj = {), (24) 

where s = tr (nij) = n r r + i^ee and we have defined k = 
Koe and A r = «oo — K rr- Also, the signal + corresponds 
to n > and — to n < 0. We can also observe from Eqs. 
(|22I24[) that under the change k — > —n, the solutions go 
asg^- j , a — )• a , ui — >• — w. 

We now discuss the magnetic flux quantization. We 
first rewrite the BPS energy density (|15|) in terms of the 
ansatz pT|) . that is, 



^ BPS = ±-S ± 

(25) 

The first term is the magnetic field whose integration 
under boundary conditions (|19I20[) gives the magnetic 
flux contribution to the BPS energy. The remaining 
three terms are total derivatives whose integration un- 
der boundary conditions (|19l20p gives null contribution 
to the total BPS energy. Thus, 



E 



BPS 



± d 2 t £ bps = ± / d 2 t B (t) 



(26) 



dt t I -— I = ±27m 



This shows that the BPS vortex solutions present energy 
or magnetic flux quantization. Next, by using the BPS 
equations (|2"2"|) and the Gauss's law (f2~4")l , the BPS energy 
density (j2"5"j) takes the suitable form 



£bps = {l-s)B 2 + + 2g 2 Q 2 + (1 + A r ) {Co') 2 , 



which is positive-definite for 

s < 1 , \ r > — 1. 

A. Asymptotic behavior 



(27) 
(28) 



The asymptotic behavior for t — > is obtained solving 
Eqs. (|22j|24p using power-series method. Thus, we attain 



9 (t) = Gt n + 



a (t) = n 



1 (1 + A r ) 



2 (1 - s) (1 + Ar) + k 2 



(t) = OJQ + 



(1 - s) (1 + A r ) + K 2 



t + 



(29) 
(30) 

(31) 
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From Eqs. (f30|) and (fT8|) . the magnetic field in the origin 
(t = 0) is given by 



B(0) 



(1 + A r ) 



(l-s)(l + A r ) + K 2 ' 
while Eq. ([3"Tj) yields the electric field ui' at i = 0, 

w'(0) = 



(1 - s) (1 + A r ) + k 2 



(32) 



(33) 



which establishes the second boundary condition for the 
field u)(t). We should note that the denominator in 
the last two equations is positive-definite due to the 
energy positivity conditions established in Eq. 



(1 — s) (1 + A. r ) + k > 0. Hence, the physical fields are 



well defined in the origin whenever conditions ([28} arc 
satisfied. 

In the sequel we study the asymptotic behavior for 
t — > +oc, for which it holds g = 1 — Sgi , a = ±5ai, 
Hi = ±<5wi, with Sgi, 5a\, Slui being small corrections to 
be computed. After replacing such forms in Eqs. (|22M24[) . 
and solving the linearized set of differential equations , 
we obtain 

6 9l ~ t-^e-P* ~ Su u Sax ~ t l ' 2 e~^, (34) 
where j3 is given as 



,3 = 



1 



2 + A r - s ± yj (A r + sf 
(1 - s) (1 + A r ) + t 



-4k 2 



(35) 



Here, (+) correspond to A r + s > and (— ) to A r + s < 0, 
such that in the limit k = we get the asym ptotic be- 
havior for the BPS uncharged vortex, (3 = y / 2/(l — s), 
see Ref. (24|. We now analyze the (3— parameter. First, 
the condition 2 + X r — s > is guaranteed because of 
the energy positivity condition (|28|) . The same holds for 
the denominator (1 — s) (1 + A r ) + k 2 > 0. On the other 
hand, the term inside the square root in the numerator, 
(A r + s) — 4k 2 , is not definite-positive. Thus, we have a 
region |A r + s\ > 2k where (3 is a positive real number, 
yielding an exponentially decaying asymptotic behavior. 
On the other hand, in the region |A r + s\ < 2k the pa- 
rameter (3 becomes a complex number with positive real 
part, which implies a sinusoidal behavior modulated by 
an exponentially decay factor. 



B. Numerical solutions 

Wc now analyze the case where /3 is a real number by 
setting A r = and s = 2k, so that j3 = y/2/(l — k). 
Consequently, the only free parameter is k, whose values 
assuring a positive-definite energy density are k < 1/2— 
in accordance with the condition ([28)) . 

In Figs. [THS]we present some profiles (for the winding 
number n = 1) generated by numerical integration of the 
equations ([22H24p using the Maple 13 libraries for solving 



the coupled nonlinear differential equations. In all figures 
the value k = reproduces the profile of the Maxwell- 
Higgs vortex [2| which is depicted by a solid black line. 
The legends given in Fig. [T]hold for all figures. 




FIG. 1: Scalar field g(t) (Solid black line, k 
solution for the Maxwell-Higgs model). 



0, is the BPS 



Figs. [T] and [5] depict the numerical results obtained for 
the Higgs field and vector potential, whose profiles are 
drawn around the ones corresponding to the Maxwell- 
Higgs model. These profiles become wider and wider for 
k < and increasing \n\, reaching more slowly the re- 
spective saturation region. Otherwise, for < k < 1/2 it 
continuously shrinks approaching the minimum thickness 
when k approaches to 1/2. Moreover, the vector poten- 
tial profile displays a novelty: for < k < 1/2 it assumes 
negative values over a small region of the radial axis (see 
insertion in Fig. [2]), which is obviously associated with 
a localized magnetic flux inversion. This inversion be- 
comes more pronounced for k values near to 1/2. Also, 
the region presenting localized magnetic flux inversion is 
a little shifted to the origin when k tends to 1/2. 




FIG. 2: Vector potential a(t). 

Fig. [3] depicts the magnetic field behavior. The pro- 
files are lumps centered at the origin whose amplitudes 
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are proportional to (1 — k)~ 2 . For n < and increas- 
ing values of |«|, the magnetic field profile becomes in- 
creasingly wider with continuously diminishing ampli- 
tude. For < K < 1/2, the profile becomes narrower and 
higher for an increasing k, reaching its maximum value 
for K = 1/2. A close zoom on the profiles corresponding 
to k closer to 1/2 (see insertion in Fig. [3]) reveals that the 
magnetic field flips its signal, showing explicitly localized 
magnetic flux reversion. 




1 2 t 3 4 5 

FIG. 3: Magnetic field B(t). 

In a first view, Fig. [3] seems to show that the mag- 
netic flux varies with the value of k. Note, however, that 
the magnetic flux is calculated by integrating 2irtB(t), 
as showed in Eq. (|26|) . An explicit plot of the func- 
tion 2irtB(t) clearly shows that for every k < 1/2 the 
magnetic flux is the same independently of the localized 
magnetic field reversion. This result is in accordance with 
the properties of the BPS-vortex solutions. 



CO' 1 




FIG. 4: Electric field u'(t). 

Fig. [4] shows that electric field profiles also are lumps 
centered at the origin with amplitudes proportional to 
k/(1 — ft) 2 , having a minimum value for k — — 1 and max- 
imal value for k = 1/2. As it occurs with the magnetic 



field, the electric profiles become narrower and higher 
while « increases tending to 1/2. Now, the difference is 
that, for ft < 0, the profiles become negative (as predicted 
after BPS equations). A close zoom along the t— axis (see 
insertion in Fig. [4]), however, reveals that the electric 
field undergoes inversion both for positive and negative 
values of k. Such reversion is ubiquitous in all profiles. 




FIG. 5: Energy density £{i). 

The BPS energy density (f27|) is exhibited in Fig. [5j 
having profiles very similar to the magnetic field ones, 
being more localized and possessing a higher amplitude, 
however. As the BPS energy density is positive-definite, 
no inversion regions are observed, naturally. 

By looking at the profiles of the BPS solutions it is 
observed that the spatial thickness is controlled by the 
Lorentz- violating parameter, allowing to obtain compact- 
like defects as in the uncharged case analyzed in Ref . 24 1 . 



The present model is being regarded as an effective elec- 
trodynamics, which subjected to the usual vortex ansatz, 
provides vortex solutions in a dielectric continuum [28j |. 
as already mentioned in Ref. [2~H . This interpretation al- 
lows to consider Lorentz- violating parameters with mag- 
nitude above the values stated by the known vacuum 
upper-bounds. 

These charged vortex configurations arc endowed with 
several interesting features, as space localization (ex- 
ponentially decaying behavior), integer magnetic flux 
quantization, magnetic flux and electric field reversion. 
Specifically, the localized magnetic flux inversion is a 
very interesting phenomenon, with sensitive appeal in 
condensed matter superconducting systems. Recently, 
a magnetic inversion was reported in the context of frac- 
tional vortices in superconductors described the two- 
component Ginzburg-Landau (TCGL) model [29]. In it, 
the magnetic flux is fractionally quantized and delocal- 
ized, presenting a 1/r 4 decaying behavior, and a subtle 
reversion. Such scenario, however, differs from the one 
described by our theoretical model, which provides a lo- 
calized magnetic flux of controllable extent, exhibiting 
flux inversion just for the parameters that yields nar- 
rower (compactlike) profiles. 
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One should still note that the set of parameters here 
analyzed is only one of the many theoretical possibilities 
available. So, it can exist a set of parameters for which 
the magnetic flux might undergo a more accentuated re- 
version, reaching more appreciable flipping magnitudes. 
Such behavior could be associated with a complex (3 pa- 
rameter, yielding oscillating solutions which become less 
localized (and more similar with the ones of Ref. [29|). 
Another interesting way is to investigate this theoreti- 
cal system in the presence of Lorentz- violating terms in 



the Higgs sector. Such analysis are under progress with 
results being reported elsewhere. 
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